Clustering dynamics of Lagrangian tracers in free-surface flows 
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We study the formation of clusters of passive Lagrangian tracers in a non-smooth turbulent 
flow in a flat free-slip surface as a model for particle dynamics on free surfaces. Single particle and 
pair dispersion show different behavior for short and large times: on short times particles cluster 
exponentially rapidly until patches of the size of the divergence correlation length are depleted; on 
larger times the pair dispersion is dominated by almost ballistic hopping between clusters. We also 
find that the distribution of particle density is close to algebraic and can trace this back to the 
exponential distribution of the divergence field of the surface flow. 
PACS: 83.50.Lh, 47.27.Eq, 47.27.Qb 
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The Lagrangian evolution of passive tracers in turbu- 
lent flows has attracted considerable attention recently. 
Particle tracking techniques have been developed that 
allow for a detailed observation of their motion even in 
high- Reynolds number turbulence . The statistics of 
few particle clusters has been used to obtain informa- 
tion on the Lagrangian statistics of the flow field [Q. 
Simplified models of passive scalars evolution, e.g. the 
Kraichnan model with its delta-correlated random veloc- 
ity fields, have provided important insights into the origin 
of intermittency corrections to scaling laws J^] . The ad- 
vection of particles that are not neutrally buoyant gives 
rise to clustering and this process has been suggested to 
be essential for the formation of rain jfj). Similar clus- 
tering phenomena should appear for bubbles or inertial 
particles in turbulent flows |7 ] . The problem we consider 
here is the advection of Lagrangian particles in a fiat 
free surface above a turbulent volume flow. Previous ap- 
proaches to the particle advection in such (compressible) 
flows include random maps S and Kraichnan models 
with prescribed smooth || and non-smooth [jl0| spatial 
variations. Realistic flows have some finite time correla- 
tions, but, as we will demonstrate here, they show further 
differences: the distribution of values for the divergence 
of the flow field is exponential (and not Gaussian ||), 
the distribution of density is algebraic (and not lognor- 
mal H ) and the two-particle dispersion shows an almost 
ballistic regime for large separations. Furthermore, this 
behavior is outside the range of Kraichnan type models 
since the ratio between divergence and velocity gradient 
fluctuations is such that the surface flow belongs to a 
marginal situation where the Kraichnan models predict 
neither clustering nor exponential separation of particles 
I0§. 

The experimental realization, the dynamics and the 
properties of the flat free surface flows that we consider 
have been discussed in detail before [fil 12|. What is 
needed here is the presence of a r 2 / 3 -scaling in the in- 
ertial subrange due to connection to bulk turbulence in 
the volume below, i.e. a non-smooth flow with finite time 
correlations which goes beyond all previous analytical 



approaches. On the numerical side we integrate parti- 
cle trajectories using a bicubic spline interpolation which 
was checked by comparison with analytical examples |l5| l 
and with a direct spectral evaluation of the velocity be- 
tween the grid meshs pjjfl . The pseudospectral simula- 
tions of the Navier-Stokcs flow are based on grids with 
256 x 256 x 65 nodes, with a forcing that maintains a fixed 
energy injection rate e Jig ]. A Taylor- Reynolds number 
R\ = u 2 rrns / [v{d x u x ) r mJ i — 145 was achieved and the 
Kolmogorov length r\ = (i/ 3 /e) 1 / 4 is 0.8 grid spacings. 
The Kolmogorov time is t v = {v/e) 1 / 2 . 

A typical particle distribution that emerges from a uni- 
form initial distribution is shown in Fig. [I] (upper panel). 
The particle dynamics shows two time regimes, a quick 
clustering into elongated structures, followed by a slower 
exchange of particles between structures. Superimposed 
on the particle distribution we show the surface flow that 
can be considered as a superposition of an irrotational 
and a gradient part by the Hclmholtz decomposition the- 
orem, 



v = v. 



+ Vp = Vxq)(x,y) e z + Vip(x,y) , 



(1) 



with scalar potentials <fi(x, y) and ?p(x, y). It suggests that 
the aggregate of particles is dominated by the gradient 
field, with the particles clustering in the minima of the 
potential ip, similar to compressible Kraichnan flows [^). 
When the velocity field is projected onto the solenoidal 
part v s the particle distribution remains essentially uni- 
form and there is no clustering (lower panel of Fig. 1). 

The dynamics in this initial period is dominated by the 
exponential contraction in v p . This follows, e.g. from the 
advection-diffusion equation for a smooth density p with 
diffusivity D, 



d t p = -(V • v)p - (v ■ V)p + DAp, 



(2) 



where the divergence patches cause an exponential vari- 
ation that typically is faster than the variations of 
the other terms [[L2|. The natural time scale r = 
((V-v) 2 ) -1 / 2 w 3.57",; is approximately the lifetime of 
a divergence patch which comes out to be about 5t^ 
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(half width half maximum of the temporal divergence 
correlation function). As expected, the divergence- free 
advection contributes little to the clustering. For the dis- 
crete particles we measure the density by coarse grain- 
ing, i.e. dividing the plane into 256 x 256 grid cells and 
count the particles inside the cells. The maximal num- 
ber n max /N with total particle number, N, increases ini- 
tially exponentially, as demonstrated in Fig. g. For t > r 
the exponentially fast formation crosses over to a slower 
regime. The accumulation of more particles into larger 
clusters continues: the inset of Fig. || shows the integrated 
probability to find cells with no fewer than No particles, 
s(t,No) = J2n=N P(ti n )i f° r different values of N with 
1 < No < N. The probabilities continue to vary even 
when the maxima remain essentially constant. 

The divergence also determines the particle distribu- 
tion p{n), as shown in Fig. |3|. In Kraichnan type models 
this distribution comes out to be lognormal for smooth 
flows and delta- like for non-smooth cases ]l0| . Al- 
though we were limited to moderate resolutions and par- 
ticle numbers, our data in Fig. || are closer to an alge- 
braic distribution which is consistent with the exponen- 
tial distribution of the divergence (see Fig. |]). Neglect- 
ing the diffusion term in (0), the density increases ex- 
ponentially in the Lagrangian frame, 7(f) = p(t)/p(0) = 
exp(- f Q (V-v)dt'). For t < 1 we can simplify the ex- 
ponent to Xt where A = V • v (in units of t -1 ). If we 
assume that the density variations are faster than the 
changes in the velocity field, we have local fluctuations 
in the divergence that give rise to locally varying density 
fluctuations. If P(X) is the probability density function 
(PDF) for the divergence, then the PDF for 7 becomes 

P( 7 ) = J dX 5( 7 - c At )P(A) = ip(ln 7 /t) . (3) 

Hence, if the divergence fluctuations are Gaussian, as in 
smooth Kraichnan flows |J , the PDF of the density fluc- 
tuations is lognormal. However, in the case of the surface 
flow the divergence fluctuations have a filamentary small- 
scale structure Q (refered to as a shocklet (negative 
divergence) in compressible supersonic turbulence [p^[). 
The small scale structures appear in the PDF as expo- 
nential tails (see Fig. |). If we let P(A) = 1/ (2s)c^ I A l / s , 
where s rs 0.95 from Fig. 4, then 

P( 7 ) „ | 7 |-i-V(*t) . (4 ) 

As shown in the inset of Fig. ||, the slope —a in the tails 
of the distribution decreases with time, from about —3.0 
at t = 0.5 to —2.3 at t = 1, in good agreement with the 
prediction from (||). For much longer times the discret- 
ness of the particles shows up and the distribution ceases 
to change, for t = 21 we got a = 1.5. 

The change in behavior for times larger than about 3 
(in units of the divergence time) is connected with the 
discretness of the particles. The exponential contraction 
near the minima in the potential leads to a depletion of 



particles in the neighborhood, so that the density cannot 
increase further by accumulation once all particles that 
initially were in a region with negative divergence are col- 
lected in a cell. Estimates of the size of the cells can be 
based on the spatial correlations of the divergence field, 
C{n) = (V • v(r i0 +r,)V • v(r i0 )) with n = x, y. The typ- 
ical extension of the patches, identified from the first zero 
of the correlations, is about 20rj (inset of Fig. |]). Based 
on this decorrelation length l& w 20ry the typical maxi- 
mal number of particles is about n/N ~ I 2 ,/ L 2 ~ 0.01, in 
good agreement with the data in Fig. || 

As a second set of characteristics we consider the 
single-particle dispersion, cr(i), and the two-particle pair 
(or relative) dispersion, d(t). The first one is defined as 
the root mean square of the absolute particle displace- 
ment, a(t) = ([x(t; x , 0)— x ] 2 }]/ 2 , where (-)c denotes an 
average over the single Lagrangian particles. The second 
uses the difference Ri2(i) = xi(t;xi, o ,0) — x 2 (t;x 2 ,o,0) 
between Lagrangian particle tracks that start at x^o and 
X2 ! o and is defined as the root mean square value for all 

particles pairs, d(t) = ([Huit) — Ri2(0)] 2 )^/ 2 . In order to 
fix the dependence on initial separation we take it to be 
about 1?7, with random orientation in space. The single- 
particle dispersion reflects the influence of flow structures 
at different scales on the particle motion and the rel- 
ative dispersion can detect the clustering, an interpar- 
ticle property. In two-dimensional incompressible flows, 
the limiting cases for both quantities are well-known [fl7| . 
Both quantities have a ballistic regime, ~ t, for short 
times when particle distances lie within the viscous sub- 
range. For times much larger than the Lagrangian inte- 
gral time scale, Tl, correlations can be expected to have 
decayed, and the relative or single particle motion be- 
comes statistically independent and both dispersions in- 
crease diffusively as in an uncorrelated Brownian motion, 
i.e. ~ t 1 / 2 . For intermediate times, anomalous scaling, 
~ t 13 with 1/2 < j3 < 1 has been observed. For an inverse 
Kolmogorov cascade, pair dispersion scaling exponents 
were found to be close to the classical Richardson value 
of 3/2 [jlSj in numerical simulations (l!]j2(J as well as 
experiments plf , but, e.g., sensitive to initial pair sepa- 
ration. In the case of single particle dispersion transient 
trapping of tracers in coherent vortex structures [fl7| af- 
fects the value of (3. 

Results on a(t) and d(t) for the full surface flow and 
the solenoidal part alone are shown in Fig. ||. The integral 
length scales are Tl/t^ — 9.1 for v and slightly shorter, 
Tl/tt) — 7-8, for v s . In all cases we do observe the initial 
ballistic regime up to T/,. The single-particle dispersion 
crosses over to the Brownian regime, a(t) ~ t 1 / 2 , for 
t > Tl, in both flow fields as indicated in the upper panel 
of Fig. [|. We connect this behavior to the fastly varying 
divergence patches |l2j that cause a kind of stochastic 
sweeping of the tracers. 

For intermediate times pairs separate superdiffusively 
like d(t) ~ t@ with an exponent of about 1.6, a value that 
is close to the Richardson prediction, rffl(t) ~ t 3 / 2 . Small 
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differences may aiso be attributed to the additional fact 
that the surface flow was found to have larger intermit- 
tency corrections than the associated three-dimensional 
volume turbulence jnj. Similar scaling behavior is ob- 
served for the pair advection in the solenoidal part only 
(see the lower panel of Fig. || for both) . 

While this anomalous scaling continues for advection 
in v s to even larger times, we find a change of the pair 
dispersion to an almost ballistic behavior of d(t) ~ t 9 
for t > 50t v , that does not seem to cross over into a 
Brownian regime. Microscopically, this means that while 
one particle follows its partner within a pair, pair corre- 
lations decay more slowly and d(t) grows more rapidly 
than in the Brownian case. Such almost ballistic scal- 
ing was also found for the single-particle dispersion in 
the strongly compressible one-dimensional Kuramoto- 
Sivashinsky equation p2| . The difference to our case 
might be caused due to the dimensionality of the problem 
and the specific character of our surface flow. Superim- 
posed on this process is a chaotic component that comes 
from v s and causes exponential separation. Another sep- 
aration mechanism is a breaking of larger particle clus- 
ters due the rapidly emerging flow, i.e. due to the rapidly 
changing divergence patterns. Clustering and separation 
were found to be competing processes that cause anoma- 
lous diffusion on a longer transient phase of the evolution. 

The exponential concentration described above has im- 
plications for the dynamics of (inertial) particles that are 
not density matched with the fluid in which they move. 
The relations by Maxey and Riley |23| for their motion 
implies that the velocity field of the particles is not diver- 
gence free. The particles will then cluster exponentially, 
as in Eq. (|3[). With a view towards the formation of rain 
there is a uniform condensation of droplets from ther- 
modynamic nucleation and then an exponential cluster- 
ing to form larger drops, which then fall to the ground 
as rain drops. The sizes of clusters thus range from the 
small scale droplets to the size of rain drops and their 
size distribution thus reflects the distribution of diver- 
gence fluctuations by Eq. (||). 
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FIG. 3. Probability density p(n) for different times during 
the cluster evolution. Straight lines indicate fits with an alge- 
braic law. The inset shows the exponents a of p(n) ~ rT a as 
a function of time t/r. The dotted line follows from (4). 
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FIG. 1. Distribution of 36000 tracers and the instantaneous 
velocity field. Upper panel: full flow v. Lower panel: solenoidal 
flow v s . The snapshots for both cases were taken at t/r v = 21 
after the start. In order to highlight the tracer patterns the 
underlying flow fields are shown in one half of the box. 
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FIG. 2. Maximum particle number per cell n max /N vs. t/r 
for particle numbers 36000, 18000, 9000, 4500, and 2250 (from 
bottom to top). Dashed line is for v s and solid lines are for 
v. The inset shows the temporal evolution of the integrated 
probability s(t,No) with values of No to the right. 



FIG. 4. PDF of the divergence of the surface flow v. The 
dotted line is an exponential fit, p(A) = 0.68 exp( — | A|/0.95) 
with A = (V ■ v)t. The inset shows the spatial and temporal 
correlation functions of the divergence field, C{x) over x/rj, 
C(y) over y/r/, and C(r) over t/t v . 



FIG. 5. Upper panel: Single-particle dispersion, a(t), as a 
function of time. Solid curve is for advection by the full surface 
flow v, while the dotted one is for advection by the solenoidal 
part v s only. The Lagrangian integral time scales for both 
fields are indicated by arrows. Lower panel: Two-particle pair 
dispersion, d(t), as a function of time. Linestyles are as above. 
The dashed horizontal lines indicate half the box size. 
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